Evidence for ?// decays into 777° and 777 
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The decays ip' 77r°, 777 and 777' are studied using data collected with the BESIII detector 
at the BEPCII e + e~ collider. Processes ip' — s> jit and ip' — > 777 are observed for the first time 
with signal significances of 4.6<r and 4.3a, respectively. The branching fractions are determined to 
be: Biip' ->■ 771- ) = (1.58 ± 0.40 ± 0.13) x 10" 6 , B(ip' -t 777) = (1.38 ± 0.48 ± 0.09) x 10~ 6 , and 
B(ip' -» 777') = (126 ± 3 ± 8) x 10" 6 , where the first errors are statistical and the second ones 
systematic. 

PACS numbers: 14.40.Gx, 12.38.Qk, 13.25.Gv 



The study of vector charmonium radiative decay to a 
neutral pseudoscalar meson P = (jr , T], rf) provides im- 
portant tests for various phenomenological mechanisms, 
such as the vector meson dominance model (VDM) [l|- 
[H, two-gluon couplings to qq states 0, mixing of r) c — 
@i HI > and final-state radiation by light quarks [l| . 
Direct contributions from the continuum through a vir- 
tual photon: e + e~ — > 7* — > jP are relevant to the decays 
of J/ip, ip' and ip" to 7P as discussed recently in Ref. Q. 
Furthermore, the possible interference between the char- 
monium decays and continuum process may play a key 
role in understanding the difference between J / ip and ip' 
decays into 7P 0]. 

For P = 77 and 77', the ratio Rj/f = B(J/tp —> 
-fi])/B(J/ijj — > 777') can be predicted by first order per- 
turbation theory [lj. The analogous ratio (R^') can be 
defined for ip' radiative decays into 77 and 77', and R^,> ss 
Rj/^j is expected Recently, the CLEO Collaboration 
reported measurements for the decays of J/ip, ip' and ip" 
to 7P [8], and no evidence for ip' —¥ 777 or 771" was found. 
Therefore they obtain R^> <C Rj/ip with R^> < 1.8% at 
the 90% C.L. and R J/ip = (21.1 ±0.9)% [|. Such a small 
R^,i is unanticipated, and it poses a significant challenge 
to our understanding of the cc bound states. 

The decay ip' — > jir° is suppressed in QED because the 
photon can only be produced from final state radiation 
off one of the quarks. It has also been described via the 
strong process ip' -> ggg -> p*7r°, p* -> 7 in the VDM 0|. 
In Ref. Q , the contribution from -0' — !• 7'* — >■ 771-° is calcu- 
lated, and B(ip' jir°) w 2.19 x 10~ 7 is obtained, which 
is compatible to the VDM contribution and does not con- 



tradict the upper limit of 5.0 x 10~ 6 [at the 90% C.L.) 
reported by the CLEO Collaboration [|. The 7* -7-77° 
vertex was shown [9j to be characterized by a form fac- 
tor F{Q 2 ), where Q 2 = —q 2 and q is the four-momentum 
of the virtual photon 7*. By using e + e~ — > e + e~ir°, 
the form factor was measured in the CLEO [l(| and 
BABAR [ll| experiments for spacelike nonasymptonic 
momentum transfer in the range \q 2 \ = 1.6 — 8.0 GeV 2 
and 4—40 GeV 2 , respectively. The e + e~ — > tp'/j* — > jtt° 
process will be very useful in testing the form factor for 
timelike photons Q 2 = — q 2 < 

In this Letter, ip' — > 771° is studied using 7T° — > 77 
decay, ip' — > 777 is measured using 77 — s- 7r + 7r~7r° and 
7/ — > 7r 7r°7T with 7r° — >• 77, and ip' — > 777' is studied 
using 77' — > 77r + 7r~ and 77' — > 7r + 7r~?/ with 77 — > 77. 
The analyses use a data sample of 156.4 pb _1 collected 
at the ip' peak with the BESIII detector operating at 
BEPCII [lj, EH • By measuring the production of mul- 
tihadronic events, the number of ip' decays is found to 
be (1.06 ±0.04) x 10 8 Q. An independent data sample 
of 42.6 pb _1 taken at y/s = 3.65 GeV is utilized to de- 
termine the potential background contribution from the 
continuum. 

BEPCII is a double-ring e + e~ collider designed to pro- 
vide e + e~ beams with a peak luminosity of 10 33 cm _2 s _1 
at a beam current of 0.93 A. The cylindrical core of 
the BESIII detector consists of a helium-based main 
drift chamber (MDC), a plastic scintillator time-of- flight 
system (TOF), and a CsI(Tl) electromagnetic calorime- 
ter (EMC), which are all enclosed in a superconducting 
solenoidal magnet providing a 1.0 T magnetic field. The 
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solenoid is supported by an octagonal flux-return yoke 
with resistive plate counter muon identifier modules in- 
terleaved with steel. The acceptance of charged particles 
and photons is 93% over 4n stereo angle, and the charged- 
particle momentum and photon energy resolutions at 1 
GeV are 0.5% and 2.5%, respectively. 

The BESlll detector is modeled with a Monte Carlo 
(MC) simulation based on GEANT4 [H, [l6[ . evtgen [l7| 
is used to generate ip' ~^ 1^ iJVi iv' events, where the 
angular distribution of the radiative photon from ip 1 de- 
cay is l+cos 2 # in the ip' frame. The decay 77 — > Tr + ir~ir° is 
generated according to the Dalitz distribution measured 
in [l8| and 77' — > 77t + 7t~ is simulated assuming it is medi- 
ated by p° —> n + ir~ , while the decays of 77 — > 7T 7r 7r and 
77' — > ir + 7r~r) are generated with phase space. ip' decays 
are simulated by the MC event generator KKMC with 
known decays modeled by the evtgen according to the 
branching fractions provided by the Particle Data Group 
(PDG) [20|, and the remaining unknown decay modes 
generated with LUNDCHARM [l?} • 

Charged tracks in BESIII are reconstructed using 
MDC hits. To optimize the momentum measurement, 
we select tracks in the polar angle range | cos6>| < 0.93 
and require that they pass within ±10 cm from the Inter- 
action Point (IP) in the beam direction and within ±1 cm 
of the beam line in the plane perpendicular to the beam. 
All the charged tracks are assumed to be pions, and parti- 
cle identification (PID) is not required, except in 77' — > 7/? 
where the dE/dx information has been used to suppress 
QED background, most of which is from e + e~ — > e + e~7. 
Either zero or two tracks with net charge zero are re- 
quired for the final tt° /ij/rj' decay products. 

Electromagnetic showers are reconstructed by clus- 
tering EMC crystal energies. The energy deposited in 
nearby TOF counters is included to improve the re- 
construction efficiency and energy resolution. Showers 
identified as photon candidates must satisfy fiducial and 
shower-quality requirements. For the ip' — > 777 and 777' 
analyses, the photon candidate showers are reconstructed 
from both the barrel and endcap of the EMC, and show- 
ers from barrel region (| cos(9| < 0.8) must have a min- 
imum energy of 25 MeV, while those in the endcaps 
(0.86 < I cos6>| < 0.92) must have at least 50 MeV. The 
showers in the angular range between the barrel and end- 
cap are poorly reconstructed and excluded from the anal- 
yses. To exclude showers from charged particles, a pho- 
ton must be separated by at least 10° from any charged 
track. The EMC cluster timing requirements are used to 
suppress electronic noise and energy deposits unrelated 
to the event. 

Events with the decay modes shown in Table U arc 
selected. Every particle in the final state must be ex- 
plicitly found, and their vertex must be consistent with 
the measured beam spot. The sum of four- momenta of all 
particles is constrained to the known ip' mass [20j and ini- 
tial e + e~ three-momentum in the lab frame. The vertex 
and full event four-momentum kinematic fits must satisfy 
Xy x < 100 and xlc < 40 , respectively. For 77 — > 37r°, 



a looser restriction of \\c < ®® * s applied to increase 
efficiency. Further selections are based on four-momenta 
from the kinematic fit. In 77/77' channels, photon pairs are 
used to reconstruct 7r° or 77 candidates if their invariant 
mass satisfies M 77 G (120, 150) MeV/c 2 or (515, 565) 
MeV/c 2 , respectively. 



TABLE I: For each decay mode, the number of signal events 
(Ns), the number of scaled continuum background events 
(Nc) in the signal region, the number of expected background 
events from ip' decays (Nu) in the signal region, and the MC 
efficiency (e) for signal are given. The error on Ns is only the 
statistical error, and the signal region is defined to be within 
±3a from the nominal and 77' masses. 



Modes (ip' -> jX) 


N S 


Nc 


N R 


e(%) 


1p' —> 77T° 


37.4 ±9.5 


63.5 


1.8 


21.4 


ip' — > -yri(n + n~ tt°) 


8.9 ±3.6 


2.2 


0.0 


21.0 


ip' —5- 7r)(7r°7r°7r ) 


3.8 ±2.3 


0.0 


1.2 


10.7 


ip' — > 7r/[7r+7r - r)(77)] 


586 ± 25 


0.0 


4.7 


27.1 


if)' — > 777' (7T + 7T~7) 


1640 ± 44 


179.3 


111.7 


41.0 



For the ip' — > 771" analysis, the primary background 
comes from the continuum process e + e~ — > 77(7), where 
the two energetic photons are distributed in the forward 
and backward regions. We require that photon candidate 
showers lie in the barrel region of the EMC to suppress 
this background. Since ir° mesons decay isotropically, 
the angular distribution of photons from 7r° decays is 
flat in the ir° helicity frame. However, continuum back- 



ground events accumulate near cos(9de 



±1, where 



decay is the angle of the decay photon in the ir° helicity 



frame [14 1. To further suppress continuum background, 
we require | cos#decay| < 0.5. Another potentially serious 
background comes from e + e~ — > 77, in which one 7 con- 
verts into an e + e~ pair in the outer part of the MDC. 
If the track finding algorithm fails to find the track, the 
two showers in the EMC are identified as isolated photons 
without associated charged tracks. To suppress this back- 
ground, the number of MDC hits, Nhits: is counted in the 
sector between the two radial lines connecting the IP and 
the two-shower positions in the EMC. To take the EMC 
spacial resolution into account, the sector is extended by 
3.5 degrees on both sides. Because of high beam related 
background level, the hits in the inner 8 MDC layers 
are not counted in Nhits- Figure [1] (a) shows the Nhits 
distribution. As shown in Fig. [T](b), background from 
continuum event 7 conversions e + e~ — > 77 accumulates 
in the low mass region. After requiring Nhits < 10, this 
background is reduced dramatically, while there is still 
an accumulation of events at the ir° mass, as shown in 
Fig. ffl(c). 

After applying the above selection criteria, the mass 
spectra of 7r°, 77, and 77' candidates are shown in Fig. [21 
An unbinned maximum likelihood (ML) fit is used for 
each analysis to determine the event yields except for 
ip' — > 777(7r 7r 7r°). The signal probability density func- 
tion (PDF) in each mode is obtained from MC simula- 
tion. The shape for the continuum background is de- 
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FIG. 1: The distribution of (a) the number of MDC hits, 
and the two-photon invariant mass distributions for 771- final 
states (b) without and (c) with the N^its < 10 requirement. 
Solid histograms are MC simulated signal for ip' — > 771°; 
dashed histograms are the luminosity-scaled continuum data, 
and points are ip' data. The solid arrow indicates the require- 
ment on Nhits- 
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FIG. 2: Mass distributions of the pseudoscalar meson candi- 
dates in ip' —>■ 7P: (a) 77T°, (b) 777(7r + 7r~7r°); (c) 777(371°); 
(d) 7?7'[7r + 7r _ r7(77)]; and (e) 7?7'(77r + 7r _ ). The crosses are 
data, the solid histograms are MC simulated signal, and the 
dashed lines are the continuum backgrounds. Fits are shown 
as solid lines, background polynomials as dotted lines. In 
(c), the arrows indicate nominal selection criteria, while the 
dashed arrows show sidebands. In (e), the dot-dot-dash line 
shows the sum of the continuum background and the expected 
background from ip' decays. 

scribed by a second order Chebychev polynomial func- 
tion, and the yield and its PDF parameters are floated 
in the fit. The fitting ranges for 7T°,77 and rf are 0.05- 
0.30 GeV/c 2 , 0.40-0.70 GeV/c 2 , and 0.85-1.05 GeV/c 2 , 
respectively. The signal yield for ip' — > ^(tt ^ ^ ) is de- 
termined directly by counting the number of events in the 
signal region, which is in 0.51-0.57 GeV/c 2 , about 3 stan- 
dard deviations from the nominal value of the 77 mass (20j 
as shown in Fig. [2] (c), while the sideband regions are de- 
fined as 0.42-0.48 GeV/c 2 and 0.60-0.66 GeV/c 2 . The 
signal yields and the efficiencies are summarized in Ta- 
ble u 

The backgrounds remaining after event selection can 
be divided into two categories. One is from ip' decays, 



which can be studied using a sample of 10 8 MC-simulated 
inclusive ip' events. The other is from non-resonant pro- 
cesses or initial state radiation to low mass resonances, 
which can be studied using the continuum data sample 
collected at a center of mass energy of 3.65 GeV. The 
expected backgrounds from ip' decays are listed in Table 
HI where the number of background events is the num- 
ber in the signal region, which is defined as within ±3tr 
from the nominal 7r°, 77 and 77' masses. For ip' — » 77r°, 
the normalized number of events from ip' — > 77r°7r° is 
1.8 in the ir° signal range. For ip' — > ^(ir ^^), there 
are 1.2 events from the decay of ip' —> 77/(77)77(371°). 
For ip' 7?7'(77r + 7r~), the main background from ip' 
decays is ip' — > p°ir° which contributes a smooth back- 
ground. The QED backgrounds for ip' — > 777' \ir + ir~ '7) 
are from e + e~ — > 7/i + /i~ and e + e~ — > r ye + e^, and 
both of them give a smooth background under the 77' 
signal peak. For ip' — > 771°, a smooth background is con- 
tributed from e + e~ —> 77(7) events. The cross section 
for e + e~ — > 7* — > jir°('jr]) has been estimated using data 
collected at y/s — 3.65 GeV, and the upper limit on the 
cross section is less than 0.14(0.68) pb at the 90% C.L. 
Since the continuum cross section is small, we neglect 
possible interference between ip' — > 771° (77) signal and 
continuum 771"° (77) . All the backgrounds are summarized 
in Table H 

The systematic uncertainties for these measurements 
are summarized in Table [TIJ The uncertainties due to 
MDC track finding and photon detection are 2% per 
charged track and 1% per low energy photon. The un- 
certainty of detecting the high energy photon is less than 
0.25% which can be neglected. The systematic errors 
from 7T° (77) reconstruction is determined to be 1% per 7T° 
(77) by using a high purity control sample of J /ip — > ir°pp 
(J/ip — > r/pp) decay. The uncertainties due to kine- 
matic fits have been estimated using the control sam- 
ples with the same event topologies as those in the sig- 
nal cases, i.e. the same number of charged tracks and 
same number of photons. The systematic uncertainties 
due to the dE/dx requirements to identify charged pions 
in the ip' — > 777' ^77r + 7r - ) and N 7 in the ip' — > 77T are 
studied by using the control samples of J/ip —¥ pir and 
J/ip — > 7*7(77), respectively, with and without applying 
these requirements. 

In tp' —> 77T , the uncertainty due to the requirement 
on the MDC hits, Nhits < 10, is studied using a sample of 
J/ip — > 777, 77 — > 77 events. The ratios of events with and 
without the requirement on the number of MDC hits are 
obtained for both data and MC simulation. Taking the 
difference of opening angle between J/ip — > 777(77) and 
ip' — > 77T into account, the difference 3% is considered 
as the systematic error for the measurement of ip' — > jtt° 
and is due to the difference in the noise in the MDC for 
data and MC simulation. 

The uncertainty due to the background shape has 
been estimated by varying the PDF shape and fitting 
range in the ML fit. For the intermediate decays, the 
77(77') branching fractions and uncertainties from the 
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TABLE II: Summary of systematic errors (%). 



Sources 


tt u 77- 


-S> 7T + 7T 7T U 


?? -» 3tt u T]' - 


-> 7T + 7T 77(77) ?/ 


—5- 7T + 7T 7 


MDC track finding 




4 




4 


4 


Photon detection 


2 


2 


6 


2 


1 


7r°(?7) reconstruction 




1 


3 


1 




4C kinematic fit 


1 


3 





3 


2 


Background shape 


4.8 


6.4 







1 


Number of ip' 


4 


4 


4 


4 


4 


Cited branching fractions 





1.2 


0.7 


1.7 


1.7 


MDC hits 


3 










Number of photons 


4 










Total 


8.3 


9.4 


7.8 


7.0 


6.4 



PDG fit |20| are used. The total relative systematic 
errors on these measurements are 8.3%, 9.4%, 7.8% , 
7.0%, and 6.4% for ip' -> 7 7r°, iP' -> 777(7r+7r-7r°), 
V>' — > ~fr](ir ir ir ), ip' — > 7?7'[7r + 7r _ 77(77)], an d V'' ~~ * 
777'(7r + 7r _ 7), respectively, as summarized in Table ILT1 

The branching fractions of ip' decays to 7 and a pseu- 
doscalar meson are listed in Table Mil Taking the com- 
mon systematic errors into account, the combined mea- 
surements for ip' — ¥ 777, 777' modes are obtained. The 
PDG [20] values are also shown in Table IIIII With con- 
sidering the background shape uncertainty, we find the 
signal significance for ip' — > -fir ^^) to be 4.6(4.3)<r, as 
determined by the ratio of the maximum likelihood value 
and the likelihood value for a fit where the signal contri- 
bution is set to zero. 



TABLE III: Branching fractions (10 -6 ) from this analysis, 
where the first errors are statistical and the second ones are 
systematic, and the comparison with the PDG values [20(. 



Mode 


BESIII 


Combined BESIII 


PDG 


Ip' — > 77T° 


1.58 ±0.40 ±0.13 


1.58 ±0.40 ±0.13 


< 5 


ip' — > 77?(7r + 7r — 7T°) 


1.78 ±0.72 ±0.17 


1.38 ±0.48 ±0.09 


< 2 


-> 7ry(7r 7r 7r°) 


1.07 ±0.65 ±0.08 






tp' — > 7?/(7r+7r - ?7) 


120 ± 5 ± 8 


126 ±3 ±8 


121 ± 8 


— > 7?7'(7r+7i-~7) 


129 ± 3 ± 8 







In summary, we have measured branching fractions 
for ip' — > 77T°, tp' — > 777 and ip' — > 777' decays. For 
the first time, we find evidence for the ip' — > jtt and 
ip' — > jrj decays with signal significances of 4.6a and 
4.3cr, respectively. The evidence for ip' — > 771- will yield 
an experimental constraint on the 7* — > 771" vertex in 



the timelike regime at \q 2 \ = m 2 ^, [6j. For the ratio 
of 77 and 77' production rates from ip' decays, we obtain 
R 4 ,, = (1.10±0.38±0.07)%, where the statistical and sys- 
tematic uncertainties from the input branching fractions 
as listed in Table IIIII have been combined in quadrature 
after accounting for common systematic errors. This ra- 
tio is the first measurement, and it is below the 90% C.L. 
upper bound determined by the CLEO Collaboration [1] . 
The corresponding 77 — 77' production ratio for the J hp 
resonance was measured to be Rj/^p = (21.1 ± 0.9)% j8|. 
R^,' is smaller than Rj/^p by an order of magnitude. 
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